Resonant Activation Phenomenon for 
Non-Markovian Potential-Fluctuation Processes 



Tomas Novotny]^ 

Institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague, Czech Republic 

Petr Chvosta0 

Department of Macromolecular Physics, Faculty of Mathematics and Physics, 
Charles University, V Holesovickdch 2, 180 00 Prague, Czech Republic 
(Dated: August 11, 2000) 

We consider a generalization of the model by Doering and Gadoua to non-Markovian potential- 
switching generated by arbitrary renewal processes. For the Markovian switching process, we extend 
the original results by Doering and Gadoua by giving a complete description of the absorption 
process. For all non-Markovian processes having the first moment of the waiting time distributions, 
we get qualitatively the same results as in the Markovian case. However, for distributions without the 
first moment, the mean first passage time curves do not exhibit the resonant activation minimum. 
We thus come to the conjecture that the generic mechanism of the resonant activation fails for 
fiuctuating processes widely deviating from Markovian. 
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The resonant activation phenomenon first reported by 
Doering and Gadoua ^ has attracted much attention. 
After the famous seminal paper various models exhibit- 
ing this phenomenon were considered and the conditions 
under which the phenomenon is present were intensively 
studied We consider a simple generalization of the 
original model by Doering and Gadoua to non-Markovian 
switching potentials generated by arbitrary renewal pro- 
cesses. That is, we study the solutions of the following 
stochastic differential equation 



-X{t) = FQ{t)+m 



(1) 



where ^(t) is a Gaussian white noise process with 
{£.{t)£.{t')) = 2S{t - t') and Q{t) is a symmetric dichoto- 
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process generated by arbitrary waiting time density f(t) 
with the distribution function F{t) ~ ^i^dTf{T). We 
study both the stationary as well as non-stationary re- 
newal processes. To make the process stationary we 
have to replace the first waiting time density by h{t) — 
^[1 — (the distribution function H{t)), whit toi — 

J(j°° dttf{t). In non-stationary cases we take h{t) = f{t). 
The diffusion described by the random process X{t) takes 
place on a linear segment a; e (0, 1) with the reflecting 
wall at a: = 1 and absorbing wall at x = 0. 

First, we present the theory of the calculation of the 
switching-averaged Green's functions. More precisely, 
the motion within the safe domain x e (0, 1) will be 
described by the state-of-potential resolved conditional 
densities Ga/3{x,y;t), defined as 
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Gaf3ix, y, t) dx = Prob{A:(t) e{x,x + dx) and Q{t) = a \ X{Q) = y and Q(0) = /?} 

I 



(2) 



with a,(3 — ±. The subscripts specify the final (a) and 
initial {(3) state of the fluctuating potential while the ar- 
guments X, y give the final (x) and initial (y) position of 
the diffusing particle. Because of the absorbing bound- 
ary at X = the total probability in the safe domain is 
not conserved. Instead, it gradually leaks out into the 
boundary. It is helpful to consider auxiliary quantities 
called boundary channels which describe the process of 
leaking of the probability out of the safe domain into 
the absorbing boundary. We may even distinguish two 
boundary channels according to the state of the fluctu- 



ating potential at the moment of the absorption event. 
In order to describe the dynamics of the channel fllling 
processes, we introduce the boundary channel occupation 
Green's functions. These quantities are given by 

TTa/siy, t) = Vvoh{X{t) e B„ I X(0) = y and Q(0) = /3} , 

(3) 

where the Bj- denote the corresponding boundary chan- 
nels. Again, the second subscript (/?) relates to the ini- 
tial state of the potential while the first (a) specifies the 
boundary channel in question. It is convenient to write 
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these Green's functions in the form of 2-by-2 matrices 
denoted by the boldface letters G{x,y;t), Il{y;t) in the 
next. 

Assume for a moment that the potential is 
static, i.e. it is fixed in one of its two states. 
Then the above matrices which will play a role 
of unperturbed quantities are diagonal with the 
form G(0)(x,y;t) = diag{Gf{x,y;t),G^°\x,y;t)), 
n(o)(y;t) = dia.g{Trf\y;t),Tr''°\y;t)). The Green's func- 
tions (x, y; t) are given simply by the Fokker-Planck 
equations for the fixed potentials supplemented by the 
boundary conditions. The boundary channel occupation 

probabilities ttj^*' (j/; t) are given by the global probability 
conservation condition 



/■I 

TT'i\y\t) = 1 - / dxG 
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{x,y;t) , 



(4) 



which may be put with the help of the appropriate 
Fokker-Planck equations in a local form reading 



x=0 



(5) 



The last but one expression gives a clear physical insight 
as the operator — ^ ± is the probability current opera- 
tor for the corresponding slope of the potential. Namely, 
the rate of the channel filling process equals the proba- 
bility current into the absorbing boundary. 



After these preparatory steps, let us focus on the con- 
struction of the full Green's functions for fluctuating po- 
tentials. Basically, our procedure consists of three steps. 
First, we shall assume an arbitrary fixed sequence of the 
potential-switching events and we shall follow the particle 
diffusion in the corresponding time dependent potential. 
Secondly, we attribute to any such evolution the prob- 
ability weight of its realization. It is during this step 
that the properties of the underlying potential-switching 
process enter the calculation. Finally, we shall perform 
the averaging over the complete set of mutually exclu- 
sive evolutions. The averaged evolution is simply given 
by the sum over all possible evolutions weighted by the 
corresponding probabilities. This is a brief outline of 
the method of construction of trajectories introduced by 
Chvosta and Reineker in where also the full formal- 
ism can be found. The result of the procedure in our case 
is 



Git) = (1 - H{t))G^"\t) + f dti(l - F(t - ti))G(°)(t - ii) • S • /i(ti)G(°)(ti) 
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(6) 



where the symbol ' • ' denotes the matrix operator mul- 
tiplication, i.e. the matrix multiplication and the inte- 
gration over the internal spatial variables, and where the 
matrix S — ( J q ) represents the switching event. The 
full evolution is represented as a sum of processes with 
zero, one, two, etc. potential switching events. 

In the case of evaluation of the boundary channel oc- 
cupations, i.e. the boundary channel part of the Green's 
function Il{y; t), the above procedure of the construction 
of trajectories yields the proper whole set of mutually ex- 
clusive paths for the evolution of ^n(y; t), not for n(y; t) 
as one may see by a closer inspection. Repeating the av- 
eraging procedure (|^) for ^n(y;t) and bearing in mind 
(PI) we derive this simple relation between the boundary 



channel occupation and the safe domain parts of Green's 
functions 



d_ 

dt 



n(zj;t) 



d_ 

dx 



G{x,y;t) 



(7) 



This formula expresses the local conservation law of 
probability for the composite diffusion and potential- 
fluctuation process. 

The convolution structure of Green's function (|^) en- 
ables to rewrite the above complicated time integrals 
structures via the Laplace transform into a simple ge- 
ometrical series which may be formally summed up to 
the infinite order giving 
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G(z) = [(1 - ff)GW](z) + [(1 - i^)G(°)](z) • S • [/iG(°)](z) 
+ [(1 - F)G^%z) ■ S • [/G(°)](z) • S • [hG^%z) + ... 

= [(1 - H)G(")](z) + [(1 - F)G(")](2) • (1 - S • [/G(o)](z))-i • S • [/iG(")](z) 

I 



(8a) 



1 d 

n(y;z) = - ^G(a;,y;2:] 
z ox 



(8b) 



a;=0 



^*/(t)G(0)(t) mean 
Equations (pf) are 



Structures like [/G(°)](z) = 
a Laplace transform of the product, 
our main result for Green's functions of the composite 
process. Solving them we get the complete information 
about the absorption process, i.e. the full description of 
the time evolution of the probability captured in the in- 
dividual boundary channels. 

In the following, we restrict ourselves mostly to a 
reduced information concerning the boundary channel 
occupations. Namely, we will consider the asymptotic 
boundary channel occupation quantities defined as 

Paisiy) = Inn Trap{y\t) lim ZTTaf){y\z) , (ABCO), 

(9) 

and also the first moments of the boundary channels oc- 
cupation densities reading 

Tafsly) = / att — hm ^ . 

Jq at 2^0 z 

(10) 



These quantities are simply related to the mean first pas- 
sage times T±{y) for respective initial conditions Q(0) = 
± by 

r±{y) = r+±(y) + T_±(y) , (MFPT). (11) 

For the Markovian case, the waiting time distribution 
fimctions are /(t) h{t) = ^e~^*, 1-F{t) = 1-H{t) = 
g-^it _ ^ g^j^y function G{t) the following identity 

holds [fG]{z) = dte-''f{t)G{t) = nG{z + n). We 
use these properties in (0) to obtain 



G{z) 



^[(z + m)- 



(12) 



Thus, for Green's function in the safe domain G(z), we 
simply get the matrix Fokker-Planck equation valid for 
the Markovian switching process (cf. Eq. (4) in [Q). 

We may proceed further to analytically evaluate 
TTapiy; z) using ( ^ ) from the knowledge of the Marko- 
vian Green's function. With the help of identities being 
satisfied by G(z) we come to the final analytic result for 
the time evolution of the boundary occupation n(y; z) 
reading 

n(2/;z) = iG^i(0,0;z)-G^(0,2/;z) , for y e (0, 1) . 

(13) 
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FIG. 1: ABCO and MFPT for waiting times distributions 
with fast decay for t oo, namely, Fq (Markovian case) (solid 
line), Fi (long dashed line), F5 (short dashed line), determin- 
istic switching (dotted line), and two (5-functions distribution 
(dot-dashed line). The t± denote the mean first passage times 
for the respective slopes of the potential. The limiting values 
of MFPT in the static limit (the average of both r's) and 
in the infinitely fast switching limit (motion in the average 
potential) are also shown by the bars. 



This expression uses the Markovian Green's function 
Gii{x,y; z) for the diffusion in (— oo, 1) with the reflect- 
ing wall at a; = 1 only. 



Instead of quoting the full rather involved result for 
Il{y;z), we only present the physically transparent ex- 
pression for ABCO with the Doering-Gadoua initial con- 
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dition y = 1 



2/i cosh k + 

( sinh fc 
a[iF\\ ; — 



/icoshfc + aPii 



with fc = ■\/2/i + F-^. One can easily see that the 
probabiUty conservation conditions for asymptotic times 



+/3 



+ P_ 



1 are satisfied. We also verified that our 



expression leads to the famous result for r of Do- 
ering and Gadoua (lOa-c). Moreover, we give here the 
analytic expression for ABCO for the 'minus' channel 
— ^-^'^^^^--'^^^ ^ a quantity analogous to that de- 
picted in Fig. 4 of Q as a result of the Monte Carlo 
simulation (for that calculation, DG used the potential 
switching between F+ — 8 and F_ — 0) 



P_(l) 



2^fc cosh fc + 2[lF sinh fc - 2^Ffc + fcF^ 
2fc(2^coshfc + FA 



(15) 



The curves of P_(l), t(1) as functions of /i for F = 8 
are plotted in Fig. [| for reference to be compared with 
other results generated by non-Markovian switching. 

Next, we present the numerical results for several non- 
Markovian switching processes generated by various re- 
newal processes governed by waiting time probability 
densities jit) normalized (except for one which cannot 
be normalized) so that the mean switching time equals 
dttjit) = i. We evaluated the ABCO of the 



i, I.e. 

'minus' 



p-+(i)+p--(i) 



and the MFPT 



. Ta(i(X) for the symmetric 
1 consid- 



channel P- (1) 

initial condition Prob{(5(0) = ±} = ^ and y 
ered also in the Markovian case. To calculate these quan- 
tities we used a numerical solution of (^) employing the 
eigenmodes expansion of the unperturbed Green's func- 
tions. It is of interest to mention that for some specific 
waiting time distributions like r„ distributions used be- 
low, it would be possible to write down and solve in an 
analytic form the equations for Green's functions analo- 
gous to (^. 

In the first set of pictures. Fig. |^, we plot the results 
for the ABCO and MFPT for the stationary processes 
generated by waiting time distributions decaying fast for 
t ^ oo. The waiting time distributions used here cover 



the r„ distributions /r„ (i) 
with n = (Markov case) 
distribution j(,(t) = 5(t — 



(n+l)' 



-(/.t)"e- 



-(n+l)Ait 



1 and 

i) corresponding to the de- 



5, the delta- function 



terministic switching process, and a two-delta-function 



One can 



distribution /25(i) = \ b{t— ^) 

see that the qualitative features of the Markovian case are 
preserved even for the considered non-Markovian switch- 
ing potentials. Namely, both the wide resonant activa- 
tion maximum in the ABCO curve as well as the resonant 
activation minimum in the MFPT curve are present in all 
cases with the shape changes attributable to the various 
variances of the used probability densities (compare Fi, 
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FIG. 2: ABCO for waiting time distributions without mo- 
ments (power laws t'" in asymptotics) . The respective curves 
describes a = | stationary (solid line) and non-stationary 
(dashed line) cases, a — ^ stationary (dotted line) and non- 
stationary (long dash-dotted line) cases, and a = | non- 
stationary case (short dash-dotted line). The results of Monte 
Carlo simulations for Qf = | stationary case are also included. 



two-delta- functions, and F5 cases). We also performed 
the numerical simulations which fully confirmed our re- 
sults. 

Further, we performed the calculations for the wait- 
ing time distributions decaying slowly like power laws 
for large t, i.e. fa{t) oc for large t with a = |, |, |. 

The exact expressions for these densities are: fa{t) — 

3 



^/2^^(^lt) 



2/^* for a — 



3 5 



and /7/2(i) 



3\/3 



2/it ^ 



These densities do not have higher order moments for 
fc > a — 1. The distributions chosen above do not have 
moments starting from the first, the second, and the 
third, respectively. The results for these distributions, 
together with their non-stationary counterparts (for the 
a = I case, only the non-stationary result exists), are 
shown in Fig. |^. Moreover, we also present there the 
results of the Monte Carlo simulations for the a = | 
stationary case. One can see that the results are qualita- 
tively the same as in previous cases for the both a = | 
cases (stat. and non-stat.) and the a = f non-stationary 
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case. On the other hand, the a = | stationary case and 
the a = I case are quahtatively different since the reso- 
nant activation minimum in the MFPT curve is absent in 
these two cases. This shows that whenever a waiting time 
density with divergent first moment is involved (the sta- 
tionary a = I case has h{t) with divergent first moment) 
the generic resonant activation behaviour is spoiled. An 
analogous behaviour in physically different context was 
found by Barkai and Flcurov 

To conclude we have presented the calculations of 



the resonant activation phenomenon for non-Markovian 
switching potentials generated by renewal processes with 
various waiting time densities. We found that the results 
are qualitatively the same for all switching processes ex- 
cept for those which are generated by the waiting time 
densities with divergent first moment. For those pro- 
cesses the resonant activation minimum in the MFPT 
curve is not present. The method used for the calcula- 
tions may be easily extended for general potential profiles 
and different boundary conditions. 
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